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Normal metal - superconductor decoupling as a source of thermal fluctuation 
noise in transition-edge sensors 
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We have studied the origin of excess noise in superconducting transition-edge sensors (TES) with several 
different detector designs. We show that most of the observed noise and complex impedance features can be 
explained by a thermal model consisting of three bodies. We suggest that one of the thermal blocks and the 
corresponding thermal fluctuation noise arises due to the high-frequency thermal decoupling of the normal 
and superconducting phase regions inside the TES film. Our results are also consistent with the prediction 
that in thin bilayer proximitized superconductors, the jump in heat capacity at the critical temperature is 
smaller than the universal BCS theory result. 

PACS numbers: 85.25.Oj, 85.25.Am, 74.25.fc, 74.40.Gh 
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I. INTRODUCTION 

A transition-edge sensor (TES) is a thin superconduct- 
ing film that can be used as a sensitive thermometer when 
voltage biased within the normal metal - superconduct- 
ing transition region and read out with superconduct- 
ing SQUID sensors 1 . TES based devices are used as 
extremely sensitive bolometers and calorimeters to de- 
tect radiation in a wide energy range from gamma-rays 
to sub-millimeter radiation 2 , and typically the thermal 
conductance to the bath is controlled by mounting the 
TES on a thin insulating SiN membrane. Although the 
performance of these detectors is already excellent, the 
most sensitive TES devices have not yet reached the the- 
oretical limits in energy resolution. This is mostly due 
to excess noise that has been shown to be present in 
many devices^—. Several candidates for the noise sources 
have been proposed, such as thermal fluctuations within 
the TES 5 , fluctuations in the Cooper-pair density^ 9 - or 
phase-slip o 10 ' 11 , but a definitive answer is still missing. 

However, before resorting to more exotic noise sources 
to explain the data, one should be sure the known noise 
mechanisms are fully understood. These are the thermal 
fluctuations of the electrical degrees of freedom, or John- 
son noise, and the thermal fluctuations of the energy de- 
grees of freedom, usually called thermal fluctuation noise 
or phonon noise - . Both of these noise mechanisms are un- 
avoidable in bolometric TES detectors; moreover, their 
magnitude depends on the details of the device in ques- 
tion. It is therefore of utmost importance to characterize 
the detector accurately, both electrically and thermally, 
before one can understand the origin of noise. 

A useful tool in characterizing TES detectors is the 
measurement of their frequency-dependent complex elec- 
trical impedance 12 . Importantly, within the transition 
this impedance depends not only on the electrical, but 
also on the thermal circuit of the device, through the 
electrothermal feedback effect^. With the help of this 
technique, a more accurate picture of the electrical and 
thermal properties of TES sensors has emerged: First 



of all, it was realized that the dependence of the de- 
tector resistance on current, and not only on temper- 
ature, is critical for the detector response, as well 2 ^. 
This was later shown to influence the Johnson noise 
directly^, and some of the excess noise could then be 
explained as non-equilibrium Johnson noise. In addi- 
tion, it has become clear that for many detectors, the 
simplest thermal circuit of one heat capacity connected 
to heat bath through one thermal conductance is not 
adequat o 7 ' 14-17 ' 19 ' 20 . A more complex thermal circuit 
then adds new components to the thermal fluctuation 
part of the noise spectrurr£ i 21 i 22 . 

Here, we present a study of the noise and complex 
impedance of several different designs of TES devices. 
Many of the devices are based on the so-called Corbino- 
geometry TES (CorTES) 6 ?-, where current spreads out 
radially from a central contact with radius n into the 
outer contact at r Q , instead of flowing linearly [Fig. 1 
(a)]. Although it is a bit more complicated to fabricate, 
it offers advantages in modelling, as the superconducting 
and normal regions separate due to the non-uniform cur- 
rent density profile. This means that we can determine 
the phase boundary radius r& from the measured resis- 
tance R, which depends logarithmically on r\> in Corbino 
geometry: 
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where Rn is the normal state resistance of the device. 
Once we know the size of the N and S phases, their ex- 
pected theoretical heat capacities can be calculated as 
well. This level of theoretical description is not possible 
in more common square shaped TES devices. 

CorTES has also shown an excess noise component, if 
compared to the simplest thermal model. This was orig- 
inally explained by the fluctuations of the N-S boundary, 
or fluctuation superconductivity noise (FSN)£. However, 
here we show that most of the " excess" noise is simply in- 
ternal thermal fluctuation noise originating from a more 
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FIG. 1. (Color online) (a) Diagram illustrating the radial cur- 
rent distribution and separation into normal (N) and super- 
conducting (S) phases in a CorTES. rt and r are the radii 
of the inner and outer superconducting contacts and rt is 
the phase boundary radius, (b) A schematic side view of a 
CorTES with an absorber. The arrows indicate the path of 
the bias current. 
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FIG. 2. (a) Basic calorimeter thermal model, (b) Three-body 
model used in this work. 



complex thermal circuit within the device. By analyzing 
the data thoroughly, we suggest that most of the "ex- 
cess" noise is generated by the thermal decoupling of the 
spatially separated superconducting and normal regions 
of the devices. Data from typical square shaped devices 
also supports this picture: Due to non-uniformities of 
real devices and/or the lateral proximity effect^, phase 
separation can still take place, although not in such a 
controllable manner as in CorTES devices. 



II. THERMAL MODELLING 

The simplest thermal model of a TES device is shown 
in Fig. H^a), consisting of a single heat capacity con- 
nected to the heat bath. The complex impedance Z 
of a device with this one-body thermal circuit always 
traces a semicircle in the complex plane as a function of 
frequency 2 . Thus, it is experimentally straightforward 
to determine whether a more complex thermal circuit is 
required by studying the shape of the Z curve. When- 
ever additional thermal blocks are added to the system, 
Z develops bulges outwards from the ideal case ( see Fig. 
|3] (a)). The shape and size of the new features in Z de- 
pend on the heat capacities Ci of the added bodies and 
the thermal conductance links gi between them, in such 
a way that the features grow in size with both Ci and 

Fig. [3] (a) shows a typical measured impedance curve 
of a CorTES device, together with a one-, two- and three- 
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FIG. 3. (Color online) (a) A comparison of one- two- and 
three-body fits to a typical impedance data, (b) An example 
of a current noise spectrum from a CorTES device, showing 
the minor difference between fitting the high frequency noise 
with ITFN or with excess Johnson (M) noise. 



body theoretical model fits. Although simpler one and 
two-body models seem to fit some device designs by other 
group o 15 ' 17 ' 18 ' 24 , for our devices a three-body model is the 
simplest model that fits the data well, as is clear from 
Fig. [3] (a). Recently, simple three-body models have also 
been used by other groups for more accurate modelling 
of their device o 14 ' 16 . 

In this work, we have therefore used the thermal model 
shown in Fig. &b), which is the simplest three-body 
model that is physically justifiable for our devices. In ad- 
dition to the TES film with heat capacity Ctes where the 
bias power Puas is dissipated, there is a hanging thermal 
body with heat capacity C\, and an intermediate heat 
capacity Ci between the TES and the heat bath. Ci 
could be associated with the supporting SiN membrane, 
and C\ could represent the absorber on top of a TES film 
[Fig. 1 (b)], but in general it could be some other part 
of the device also. Most devices studied here do not ac- 
tually have absorbers. We call this the IH (intermediate 
+ hanging) model, and use analytical expressions for the 
impedance and noise of a TES with this thermal circuit 
presented in 2 ^. Their derivation and full theoretical dis- 
cussion will be published elsewhere 2 ^. We simply note 
that with each new thermal block, a new thermal noise 
source is also introduced to the system, arising from fluc- 
tuations in energy between the bodies through the ther- 
mal link. We call the noise due to any additional block 
internal thermal fluctuation noise (ITFN) to distinguish 
it from the thermal noise between the system and the 
heat bath (phonon noise). 



A. Intermediate body 

The intermediate body is characterized by its temper- 
ature Ti and heat capacity Ci . We emphasize that we do 
not make the simplifying assumptions that Ti = Ttes or 
that T 2 = Tbath, when the TES is biased. Associated with 
the intermediate block, we thus have four different values 
of dynamic thermal conductance, two for each physical 
link connected to Ci, evaluated at each temperature end 
of the link 2 ^. The parameters Gd yn — dP/dTrES an d 
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Ttes can be calculated from a careful analysis of the 
measured I-V characteristics, as described in Ref 26 . 

In our fitting procedure we use T 2 , C 2 and 
9tes,2(Ttes) as the free parameters. The other un- 
known conductances gTEsA T 2), 92,b(T bath ) and g 2 ,b{T 2 ) 
are then calculated using the free parameters and the 
known values of Ttes an d Gdym as2£ 

r _ 9TES,2{TTEs)g2.b{T2) (r> , 
gTES,2( 1 2) + 92,b{^2) 

To simplify the model and to minimize fitting parame- 
ters, it is also assumed that the links on both sides of 
the intermediate block have the same thermal exponent 
n, i.e. g T EsA T j) = ^T£ -1 and g 2 , h {T j ) = BT?~ l . This 
is physically reasonable if all the conductances are dom- 
inated by the phononic transport properties of the SiN 
membrane. 

Later, we will show that the intermediate body fitting 
parameters are nearly constant throughout the transi- 
tion, and only change when T^ath is changed. Also, C 2 
appears to depend on both the area and thickness of the 
membrane. This behavior is thus consistent with the as- 
sumption that the intermediate body represents the SiN 
membrane phonons. Furthermore, T 2 is usually always 
slightly below Ttes, and not close to the bath tempera- 
ture. 



B. Hanging body 

In contrast to the intermediate body, the second added 
heat capacity C\ in our model is hanging [Fig. [2] (b)], 
which means that no steady state power flows through it, 
and the average values of Ttes and T\ are equal. This 
simplifies the description in comparison to the interme- 
diate body, as only one thermal conductance gTES.i is 
required. The Joule power of the TES bias current is 
all assumed to be dissipated inside Ctes- In r °al X-ray 
and 7-ray devices with thick absorbers on top of the TES 
fil m 14 i 20 ' 24 , C\ could well be the absorber, in which case 
gTES,i describes the thermal conduction within the TES 
and the absorber. In devices without an absorber (most 
samples in this work), the idea is to model electronic 
degrees of freedom, thus gTES.i models the thermal con- 
duction within the TES film. 

The hanging block is naturally still a simplification, 
in reality there can be also a thermal coupling from C\ 
directly to the bath or to C 2 . However, the effect of the 
missing link, which would have a magnitude of the order 
of gTES,2, is believed to be minor for the devices studied 
here, where gTES,i is much larger than gTES,2, as will be 
seen later. Any effects of temperature gradients within 
the TES film are also not captured in this simple model. 
Numerical estimates of the gradients within a CorTES 
device 2 ^ have shown that the isothermal simplification is 
a reasonable assumption. 

For the CorTES devices with phase separated N and S 
regions, one may at this point wonder, if the blocks Ctes 



and Ci could actually directly represent the normal and 
superconducting regions of the TES film, respectively. 
The great benefit of the CorTES geometry is that we can 
a priori calculate the values of Ctes and C\ as a func- 
tion of bias, which is not possible for the more standard 
square shaped devices. As we later show, the fitted val- 
ues of the heat capacities and their behavior as a function 
of bias point do indeed follow a trend predicted by this 
interpretation. gTES,i is thus still the electronic thermal 
conductance, but now it can be affected by the location 
of the NS boundary, for example. Understanding of how 
9tes,i should behave as a function of the bias (phase 
boundary location) is quite sketchy at the moment, and 
it needs to be studied further theoretically. 

The hanging body C\ produces internal thermal fluctu- 
ation noise (ITFN) that appears in current noise spectra 
at the high-frequency side of the effective thermal time 
constant, just like Johnson nois o 20 ' 25 , see Fig. [3](b). In 
addition, when gTES,i is large enough, the roll-off for the 
ITFN noise occurs at very high frequencies, even above 
the electrical cut-off frequency of the read-out circuit. 
Thus, the ITFN noise produced by the hanging body 
can look very similar to the Johnson noise. Quite often 
in the past, excess noise in TES devices has been quan- 
tified by the parameter M— which is a measure of 
excess Johnson noise. Thus, one has to be very careful 
about drawing conclusions based on just the noise data: 
If the thermal circuit is not adequately characterized, one 
can misinterpret ITFN noise as excess Johnson noise. In 
our devices the time constant for ITFN roll-off is indeed 
near the electric cut-off of our read-out circuit, so that fit- 
ting the spectra either with the three-body ITFN model, 
or two-body + excess Johnson noise (M parameter) give 
almost identical results, as illustrated in Fig. G2(b). The 
only difference seen is that with the ITFN model, the 
roll-off has a steeper slope, as it is a higher order roll-off 
produced by the combination of the thermal and electri- 
cal circuits. However, the point is that in many cases the 
impedance can only be fitted with the three block model 
[see Fig. [3] (a)] , which supports the choice of ITFN noise 
over other noise sources. 



III. EXPERIMENTAL SETUP AND METHODS 

The experiments were performed in a compact home- 
made dilution refrigerator with a base temperature T = 
40 mK. The TES devices were voltage biased by a shunt 
resistor (8.9 mfi) at the sample stage, and a NIST two- 
stage SQUID 28 mounted on the 1 K stage (equivalent 
input current noise density ~ 4 pA/yTIz and geometric 
input inductance of 300 nH) served as the current am- 
plifier for the readout. The SQUID was operated in the 
flux locked loop, with dedicated room temperature elec- 
tronics designed by SRON. More details on the setup can 
be found m 2 ^. 

All studied TES devices were first characterized by 
measuring their resistance vs. temperature transitions 
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with a four-probe lock-in measurement, to find the nor- 
mal state resistance and critical temperature T c . 
Next, each detector was connected to the SQUID read- 
out, and a series of current- voltage (I-V) curves was mea- 
sured at several different bath temperatures. The ther- 
mal exponent n was extracted from the Tbath dependence 
of the I-V data, as explained in Ref. 2 ?.. Then, using the 
known n, the TES temperature Ttes and dynamic ther- 
mal conductance Gd yn at any point in the transition were 
determined from the I-V curves 2 ^. We also calculated the 
transition steepness parameter axoT.iv = {T/R)dR/dT, 
which is a dimensionless measure of the sensitivity of the 
TES 1 . 

Finally, a set of noise and complex impedance measure- 
ments were performed at a desired Tbath and at various 
bias points within the transition. The TES impedance 
was extracted from the measured circuit impedance by 
dividing out the effect of the read-out circuit transfer 
function, as described in ref^S. From the TES impedance 
curves, the low and high frequency limits Zq and Z x can 
be estimated, which can then be used to calculate the im- 
portant transition parameters a = (T / R)dR/dT\i and 
P = (I/R)dR/dI\ To as 



n Z 
a = — - 



4> z + r 







(3) 



where 0=1- {T bath /T TES ) n and Rq is the TES resis- 
tance. Additionally, aror.z — (n/<fr)(Zo — Ro)/(Z + 
Ro) = [2a + (n/</>)/3]/(2 + j3) from the impedance data 
should give^i the same value as aroTjv, so that we can 
check for consistency between the impedance and I-V 
measurements. Notice how oltot depends both on a 
and /?. 

We have performed impedance measurements between 
4 Hz and 100 kHz both by the white noise method 12 , 
where a white noise excitation is used and the broad- 
band response is measured, and by the lock-in method 
where a sinusoidal excitation is used at certain frequency 
points^. Both methods can in general be used, however, 
we have seen that the sine- wave method is more reliable, 
as the total heating power generated by the excitation is 
typically less in the sine-wave method. Data with heat- 
ing problems could be identified by the consistency check 
mentioned above: With excess heating axoT,z did not 
agree with oltot, iVt indicating a heating induced shift 
in the bias point. 

The measured impedance and noise data are finally 
fitted simultaneously to the IH model equations. We em- 
phasize that the fitting is done by eye, and free fit param- 
eters are varied manually, as high-dimensional non-linear 
least-squares fitting would be demanding to implement. 
In some cases we found that the three-body equations 
still did not explain all the observed noise, even if the 
impedance fit was good. In those cases we have quanti- 
fied the remaining truly excess noise as excess Johnson 
noise using the M-parameter. Note also that although 
here we only show the complex plane plots of Z, we also 



make sure that the real and imaginary parts fit separately 
as a function of frequency. Included in the fits, but not 
shown in the plots because of their small values, are the 
Johnson noise of the known shunt resistor, and the equiv- 
alent white input noise of the SQUID. We have also used 
the equation for the (lowest order) non-equilibrium TES 
Johnson noise 13 V£ = Ak B TR (l + 2/3) in the analysis 
for all the devices. Even though Z is measured only up 
to 100 kHz, the theory curves are always calculated up 
to 2 MHz, in order to see the high-frequency differences 
between the fits, as demonstrated in section ITVC II 



IV. EXPERIMENTAL RESULTS 

In this section we present the measured data and fit re- 
sults from several different detectors. Table Q] lists some 
of the key parameters for each device studied. Most of 
them have radially spreading current distribution as in 
the fully circulary symmetric CorTES device shown in 
Fig. [TJ However, the device labelled STES is divided into 
four equal size parallel slices (see inset, Fig. [5]), and the 
detectors labelled Slice 1 and Slice 2 are individual slices 
(Fig. ITl]) . Also, data from one traditional square-shaped 
device is also presented. All measurements were done at 
a regulated bath temperature of 60 mK unless stated oth- 
erwise. The detectors were not shielded against external 
magnetic fields, nor was earth's field compensated for. 



TABLE I. Device parameters for the TESs measured in this 
work. 

TES Ti/Au [nm] T c [mK] R N [mfi] SiN [fim] 



CorTES 


40/55 


98 


200 


0.30x800x750 


STES 


40/55 


99 


220 


0.30x800x750 


Slice 1 


71/105 


126 


166 


0.30x830x730 


Slice 2 


58/83 


162 


220 


0.75x830x730 


Square 


48/70 


156 


425 


0.75x460x410 



A. CorTES 

The data we have measured on a full CorTES turned 
out to be quite difficult to fit accurately even with the 
three-body IH model, as seen from the examples shown 
in Fig. |4] The general trends are, nevertheless quite 
well reproduced, including the very strong deviations of 
Z from the simplest circular shape and all the trends 
in the noise spectra. Typical values for the most im- 
portant fit parameters were C\ ~ 0.2 pJ/K, C2 ~ 0.35 
pJ/K, Ctes ~ 0.05 pJ/K and g TES ,i ~ 20 - 30 nW/K 
throughout the transition for the bare devices without 
an absorber. Fig. 2] also shows the data after the de- 
position of a 2 fjxa thick Bi absorber on top, as shown 
schematically in Fig. [lib). In that case, all other param- 
eters stayed about the same, except for Ci, which, quite 
reasonably, doubled to ~ 0.4 pJ/K. We should note that 
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FIG. 4. (Color online) Comparison of (a) the impedance 
(measured with the noise method) and (b) the noise of a full 
CorTES with (black crosses) and without (red open circles) 
an absorber, at two different bias points. The lines are best 
fits to the data using the IH model. No M-parameter was 
used. 



the obtained C values are reasonable if compared to es- 
timates from the detector size, except for C2 which is 
surprisingly large. This could be because of the 120 nm 
thick AlOx insulator layer separating the Nb bias lines 
PJb)] in the CorTES device, creating unwanted thermal 
links and heat capacity inside, as shown in Ref^ for a 
different device. 

Very accurate determination of how the addition of 
the absorber affected the detector characteristics is un- 
fortunately complicated, because the absorber process- 
ing changed also the transition properties, by lowering 
the T c from 98 mK to 85 mK, and by broadening the 
transition. For the bias values shown in Fig. |4l oltot,iv 
decreased from ~ 300 to ~ 50 at R/Rn = 0.6, whereas 
at R/Rn = 0.1 it stayed constant at ~ 150. These 
changes influence the detector responsivity, which means 
that noise and Z are affected not only through the in- 
crease of C\ but also directly. The large difference of 
noise at R/Rn = 0.6 can be attributed mostly to the 
change in a, whereas at R/Rn = 0.1 the lowering of 



noise is consistent with the lower T c of the device after 
the addition of the absorber. 



B. STES 

We have also measured a set of data from another ab- 
sorberless CorTES device, shown in Fig. [5j It was oth- 
erwise nominally identical to the full device discussed 
above, but was cut into four parallel slices (dubbed 
STES), with an SEM image shown in the inset of Fig. 
in^b). In this case, the impedance fits to the IH model 
look nearly perfect, and the noise data agrees also quite 
well without the need of an empirical M parameter, ex- 
cept at the very highest frequencies near the cut-off (rea- 
son unknown at the moment). Achieving better fits than 
in the full CorTES case was easier, most likely because 
the STES device has a lower responsivity (lower max- 
imum a ~ 200 ). In Fig. [6] we plot the obtained fit 
parameters of interest as a function of the bias point 
R/Rn- The heat capacities, Fig. [S] (a), do not vary 
strongly, and their values are generally consistent with 
the CorTES results. C2 ~ 0.55 pJ/K is again high, even 
higher than in the CorTES device. To compare the data 
with the suggestion that Ctes and C\ are associated 
with the normal and superconducting regions, respec- 
tively, we also show theoretical curves corresponding to 
the calculated heat capacities of the normal and super- 
conducting phase regions, using the location of the phase 
boundary calculated from Eq. [T] the known film thick- 
nesses and Sommerfeld constants jTi — 330 J /K 2 m 3 and 
lAu = 65 J/K 2 m 3 from literature*! 3 ". In addition, the heat 
capacity of the S region was calculated using two different 
assumptions: (a) the jump at T c is given by simple BCS 
theory as 1.43CV (dash-dotted line), and (b) the jump 
is slightly suppressed due to proximity effect according 
to Refill. We see that C\ is clearly consistent with it 
being associated with the S region, however Ctes is a 
bit elevated compared to the simplest theory. Moreover, 
lower in the transition at R/Rn < 0.3, there seems to 
be a trend that C\ increases and Ctes decreases. This 
trend is very clear in the data, as forcing the values of 
C\ and Ctes f° r R/Rn < 0.3 to equal the values at 
R/Rn = 0.3 produces very poor fits, as shown in Figs. [5] 
and [6] 



C. Slice TES 

The usual CorTES devices have an AlOx insulator 
layer between the bias leads. Furthermore, from the high 
values of C2 obtained above, and from previous measure- 
ments in Ref^ 2 -, we suspect that the AlOx layer could be 
responsible for the high value of GV Therefore, we also 
studied devices where only a quarter "slice" of the full 
CorTES disk is retained: this way the insulator layer is 
not required, as the center contact lead can come from 
the opposite side (see inset Fig. [TTj) . However, we still re- 
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FIG. 5. (Color online) Measured (symbols) and fitted (lines, 
IH model) (a) impedance (measured using the lock-in method 
up to 100 kHz) and (b) noise in the STES device. Decreas- 
ing bias makes the impedance curve smaller and noise level 
higher. Bias values shown are R/Rm = 0.1, 0.2, 0.3, 0.4, 0.5, 
0.8, rest are omitted for clarity. Dashed lines show fits with 
alternative parameter values shown in Fig. [6] as crosses. 

tain a geometry that promotes a phase separation. Here 
we report on measurements on two such samples, with 
different Au/Ti thicknesses, leading to different T c and 
Rn values (Table H]). The Au/Ti layer thicknesses were 
increased compared to the full CorTES devices, in order 
to keep Rn approximately the same. Also, the thick- 
ness of the SiN membrane (300 nm) was chosen to be 
the same as in the usual CorTES devices for sample 1, 
whereas sample 2 had much thicker SiN (750 nm). 



1. Sample 1 

The first slice TES has relatively low values of a ~ 
60 — 80, with the result that the measured impedance 
data, shown in Fig. [Jj lacks any striking features. The 
data is therefore too easy to fit: That is, we can find 
several ways to fit the data (both impedance and noise), 
and choosing the best fit is not simple based on this data 
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FIG. 6. (Color online) Parameters obtained from the fits to 
the data in Fig. [S]as a function of the bias point. (a) The heat 
capacities, and (b) the thermal conductances. Gdyn is fixed 
by the I-V data and is not a fit parameter. Theory curves 
show the expected dependence of the normal region (dashed 
line) and the superconducting region with two different mod- 
els for the jump at T c : full jump (dash-dotted) and suppressed 
according tc— (solid line). Crosses show alternate parameter 
values used for "bad fits" shown in Fig. [S] Inset in (b) is an 
SEM image of the STES. 



alone. To illustrate this, we have fitted the same data 
in three different ways, with the resulting fitted curves 
almost identical up to 100 kHz, but with significant dif- 
ferences in the obtained parameter values. The three 
cases we consider are: Case 1: we force Ctes to the the- 
oretically calculated value and let the other parameters 
vary freely. In case 2, Ctes is not fixed but also free (as 
in all the fits for the CorTES and STES devices above). 
Finally, in case 3, we set C\ — 0, reducing the thermal 
model to two bodies: the TES and the intermediate. 

From the CorTES sample results above, we still ex- 
pect that the three-block model is required to describe 
the sample physics, thus we first discuss the comparison 
between cases 1 and 2, and only later comment on the 
case 3 fits. Fig. [Jj shows a comparison of cases 1 and 2 
with the data. It is clear that both cases fit the data well, 
and deviate from each other only at high frequencies close 
to 100 kHz, with case 1 making a more pronounced kink 
in the calculated impedance and slightly weaker roll-off 
for the noise. 

Looking at the obtained parameters plotted in Fig. [51 
we find case 1 quite interesting. With Ctes forced to 
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the theoretical values, C\ follows the calculated theoret- 
ical value of the reduced BCS heat capacity^!, support- 
ing again the proposed picture of thermal decoupling be- 
tween the N and S phases. As expected, due to the lack of 
AlOx, C2 ~ 0.1 pJ/K is much reduced compared to the 
CorTES and STES samples, even though the SiN mem- 
brane dimensions are the same. Case 2 reproduces all the 
trends, but now Ctes is higher than the theory (as in the 
STES device) and C\ lower, but with their sum approx- 
imately constant between the two cases. We also note, 
how for both cases C\ grows fast below R/Rn < 0.2, 
similar to the STES device. 

On the other hand, comparing the two cases in terms 
of grES,i fit values produces a striking difference: The 
trends are quite clear, but opposite for the two cases. 
Thus, it is hard to draw solid conclusions, except that one 
can see that the highest values obtained ~ 100 nW/K are 
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consistent with the increased TES bilayer film thickness, 
if compared to the STES results. The trend in case 1 
of decreasing <7tes,i with decreasing R is of course con- 
sistent with the fact that thermal conductance starts to 
decrease in the superconducting phase, and that the size 
of the S phase increases when going down in bias. 

In cases 1 and 2 we did not need to use the M- 
parameter (excess noise) in the upper part of the tran- 
sition, but had to include it in the lower part, in case 
1 below 30 % bias and in case 2 already at 50 %. This 
is in contrast to the CorTEs and STES fits, where M- 
parameter was not used at all (those fits were not as 
accurate, though). The fitted values of M are shown in 
Fig. [Tol a). For this sample, we also found that in a 
narrow bias range Rj Rn = 0.3 to 0.35 there was an en- 
hancement of the mid- frequency (1 kHz - 30 kHz) noise 
that could be explained neither by the thermal circuit 
nor an M- noise component. Curiously, this bias range 
corresponds to the sudden drop in the a parameters, as 
can be seen in Fig. EJc)). 

To study Case 3, we show in Fig. [5] representative 
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impedance fits at two bias points (solid lines), with case 
1 fits included for comparison. Again, if we only look at 
frequencies up to the measurement maximum, we do not 
see any difference in the fits. The difference in the noise 
fits (not shown) is also limited to a small variation in the 
roll-off. In case 3 we naturally do not have the high fre- 
quency ITFN component present because C\ is missing, 
and thus a large part of the noise at high freqencies has 
to be accounted for by excess noise with M parameter, 
at all bias points. The effect of the analysis case on the 
obtained Af-parameter values vs the bias point is shown 
in Fig. [TU] (a). In addition, the chosen fitting case has a 
big effect also on (3, as it is proportional to [Fig. [10] 
(b)]. 

The plots remind us that many parameters strongly 
depend on the chosen thermal model and how it is inter- 
preted. If we do not know with certainty which of the 
possible thermal models is correct, we should not jump 
to conclusions about the nature of the excess noise based 
on the dependence of M parameter on bias, magnetic 
field etc. In other words, the Z data may look like it fits 
a simpler one- or two-block model, but the underlying 
true model could still be different. For example, if one is 
developing a detector where /3 is an important parame- 
ter, mistakes can be made in TES design if the decisions 
are based on /3 from a wrong fit. For the case of this 
particular device, we believe that Cases 1 and 2 are more 
accurate because of the evidence from the other samples 
studied here, where two-block model fits are impossible 
to the impedance data. 

Another interesting point to notice is that even though 
we have three different sets of values for Ctes an d C\, 
the sums Ctes + Ci are essentially constant between the 
different fitting cases, as shown in Fig. fTUTc). The values 
are also consistent with the modified BCS theory^!, ex- 
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cept at the lowest bias points, where the rapid increase 
mentioned before is clearly seen. 

At this point we wish to speculate about the observa- 
tions at low bias points. A possible scenario for the ef- 
fects seen (increasing C, increasing Ai-noise) could well 
be related to vortex physics. If vortices are not pinned, 
they could contribute to heat capacity and also gener- 
ate excess noise by the so called phase-slip shot noise 
mechanism3ii. The sudden drop in cxtot,iv is also cor- 
related with the onset of the need for Af- parameter in 
Case 1. The lower a would then correspond to extra re- 
sistance in the transition, caused by the vortex motion. 
In this light we can perhaps also understand why the full 
CorTES devices do not require M to fit their noise, as 
the extra Nb top layer over the TES film may pin vor- 
tices more strongly, and prevent the generation of excess 
phase slip noise. Another option is that some of the ex- 
cess M-noise could be generated by the FSN mechanism 8 . 
However, FSN theory in its current state cannot explain 
the increased heat capacity or lowered a. 



2. Sample 2 

The geometrical design of slice TES sample 2 was iden- 
tical to sample 1, however the SiN and TES layer thick- 
ness are different (Table [TJ) , resulting also to a higher T c . 
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FIG. 11. (Color online) Measured data (symbols) and fits 
(lines) for slice TES 2 at T bath = 60 mK. Insets: (a) SEM 
image of the TES, (b) Total a calculated from I-V. 



In contrast to sample 1, the second slice TES featured 
a more complicated transition with very large peaks in 
ctTOTjv as shown in the inset of Fig. [TTT b). The high 
values and large variation in oltot.iv produces compli- 
cated impedance features and high noise levels. Sample 

2 was also measured at four different bath temperatures: 
140 mK, 110 mK, 85 mK and 60 mK. Fig. HU shows the 
data for a few representative bias points measured at a 
bath temperature of 60 mK. The fits were done with the 
full IH model using Case 2, keeping Ctes free. Now Case 

3 is again out of the question due to the complex shapes 
of the Z-data. Case 1 was tried but did not produce as 
good fits as Case 2. 

Again, the main features of both Z and noise data are 
captured with the IH model. However, at bias points 
where a is very large (for example R/Rn — 0.47) it was 
more difficult to find good fits, and determination of Z x 
is more uncertain. The inevitable result of this is more 
scatter in the fitted parameters. In Fig. [T2]we have plot- 
ted all the results from all four bath temperature runs in 
one graph, as Ctes, Ci and gTES,i should not depend 
on Tbath at all according to the thermal model (C2 de- 
pends on T ath through its influence on T2, thus we show 
only Tbath = 60 mK fits). Comparing to Sample 1, we see 



that Ctes and C\ have the same trends, and only slightly 
higher values, consistent with the increase of T c . On the 
other hand, C 2 is three times larger. This again supports 
the picture that C2 originates from the SiN membrane, as 
it is 2.5 times thicker for Sample 2. gTES,\ is of the same 
order of magnitude, and does not follow a monotonous 
trend throughout the whole transition. Most of the noise 
data did not require any additional M-noise, all the fits 
shown in Fig. [IT] are without the M parameter. This, we 
speculate, could be because of the much higher a com- 
pared to Sample 1: The themal noise grows so much that 
is swamps any possible additional noise sources. 

The data for gTES,2 was not shown because it depends 
on the bath temperature. For comparison with sample 1, 
at 60 mK sample 2 had cites, 2{Ttes) = 1<1 nW/K and 
G dyn = 0.97 nW/K at 0.5 R/R N . 

As can be seen for bias point R/Rn = 0.2 in Fig. [TT] 
the measured noise sometimes develops a peak near the 
high-frequency roll-off, when the TES is biased low in the 
transition. This indicates electrothermal oscillations and 
it is predicted by our thermal model, as shown by the 
fit. It arises due to a resonant-like interaction between 
the electrical and thermal circuits: When C\ starts to 
decouple from a TES that has a very small heat capacity 
and large enough a, the TES response becomes oscil- 
latory. The increased responsivity means that all noise 
sources develop a peak, including the Johnson noise, for 
example. 
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D. Square TES 

Finally, we show also results for a simple square TES. 
The measured detector was a bare 300/im x 300/im TES 
with Nb bias lines and no extra features added (more 
sample parameter details in Table IJ). It has a smoothly 
changing axoTjv with a fairly low maximum value of 50. 
In Fig. [T3] we plot the measured and fitted impedance 
and noise data. Just as for all the other geometries, our 
square TES clearly has a large high frequency noise com- 
ponent. Again, the shapes of the impedance curves shows 
that the model without a hanging block (case 3 discussed 
in section lTV C II) will not work here either. The IH model 
again produces a very good fit for both Z and noise, 
which gives us confidence that the noise in the square 
TES can be explained by the same mechanism as in the 
other geometries. Note that no M parameter was needed 
in these fits. 

Fig. [14] shows the relevant parameters from the fits. 
We again observe a decreasing Ctes and an increas- 
ing C\ , consistent with our model of thermal decoupling 
of the N and S regions, and comparable values to the 
Corbino-like samples. To calculate some estimate for the 
N and S phase heat capacities, we take a crude model of a 
linear transition, where the ratio of the volumes is given 
directly by R/Rn, so for example in the middle of the 
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FIG. 14. (Color online) Parameters from the fits in Fig. 1131 
Lines show the simple estimates for Ctes (dashed) and Ci 
(solid) discussed in the text. 



transition we would have exactly half of the device in the 
superconducting state. This model produces a correct 
order of magnitude for Ctes an d Ci, but naturally no 
real agreement, although the parameter values do seem 
to change linearly with R, as our overly simple model 
predicts. The values of thermal conductances are again 
consistent with the Corbino devices, and <?tes,i does not 
show any simple dependence on R. 



V. CONCLUSIONS AND DISCUSSION 

We have fabricated several different geometrical de- 
signs of TES devices, and studied their impedance and 
noise properties. All data are consistent with the picture 
that a three block thermal circuit is required to accu- 
rately describe the data. In all devices, the fitted heat 
capacity values follow trends which suggest that one of 
the extra heat capacity blocks arises from the insulating 
films such as the SiN membrane and, if present, the AlOx 
layer used in some devices. The other two blocks seem 
to represent the normal and superconducting regions of 
the TES film itself, with the observation that the finite 
value of thermal conductance within the TES film (or 
a possible thermal boundary resistance between the N 
and S phases), in combination with the large value of 
the heat capacity in the superconducting phase near T c , 
leads to fairly large observable features in the impedance 



11 



and noise data at high frequencies. In most cases, this 
internal thermal fluctuation noise component explains all 
of the "excess noise" in the devices. However, in some 
devices indications of additional noise sources remained, 
especially low in the transition. 

In our Ti/Au devices the resistivity of the TES film 
is such that grES,i is "low enough". By this we mean 
that the resulting cut-off frequency of the ITFN noise 
due to the hanging block is near the electrical cut-off 
of the readout, so that we can observe its effect in the 
measured data. Our data, therefore, agrees with the pre- 
vious suggestion that a large ITFN noise exists in Ti/Au 

In contrast, the lower resistivity Mo/Au or Mo/Cu de- 
vices usually have a lower excess noise level that is white 
voltage noise, in other words noise that looks like ex- 
cess Johnson noise^^. Our model could perhaps explain 
part of this difference between Ti based and Mo based de- 
vices through the higher thermal conductance grES.i of 
the Mo/Au and Mo/Cu TES films. If g T ES,i grows, the 
ITFN noise level will fall, and the cut-off is pushed to a 
higher frequency. Thus if the ITFN roll-off happens much 
after the readout cut-off, ITFN noise will look like excess 
Johnson noise in the frequency range of the measurement. 
More measurements on Mo devices are required to test 
this hypothesis. 

The results reported here indicate that in an optimal 
TES design, the internal thermal conductances within 
the metallic parts of the device (TES film, absorber) 
should be maximized. As the actual phase separation 
between the superconducting and normal regions may 
be impossible to prevent, the effect of the superconduct- 
ing heat capacity is minimized if the heat capacity of 
the normal part of the TES is larger. We note that the 
leading designs in terms of reported energy resolutions^ 
employ added normal metal features, and thus we spec- 
ulate that the reduction of their noise could originate, 
at least partly, due to the increased normal region heat 
capacity. 
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